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Abstract 

We study a physical model for the interaction between general inclusions bound to fluid mem- 
branes that possess finite tension 7, as well as the usual bending rigidity k. We are motivated by an 
interest in proteins bound to cell membranes that apply forces to these membranes, due to either 
entropic or direct chemical interactions. We find an exact analytic solution for the repulsive interac- 
tion between two similar circularly symmetric inclusions. This repulsion extends over length scales 
~ yj K/7 and contrasts with the membrane-mediated contact attraction for similar inclusions on 
tensionless membranes. For non circularly symmetric inclusions we study the small, algebraically 
long-ranged, attractive contribution to the force that arises. We discuss the relevance of our results 
to biological phenomena, such as the budding of caveolae from cell membranes and the striations 
that are observed on their coats. These, and other, "gnarly buds" may prove fascinating to study 
further. 

PACS numbers: 87.16.-b Subcellular structure and processes 87.16.Dg Membranes, bilayers, and vesicles 
82.70.Uv Surfactants, micellar solutions, vesicles, lamellae, amphiphilic systems 
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I. INTRODUCTION 



A significant proportion of all proteins in a typical eukaryotic cell are membrane proteins. 
These are found anchored to cell membranes. Many of these carry out tasks such as signal 
transduction, pore or channel formation, cytoskeletal binding etc[l]. Others are involved 
in endocytosis and exocytosis. In particular it is now understood that the formation of 
clathrin coated pits is driven by the controlled geometric aggregation of clathrin which exert 
corresponding forces on the cell membrane [2, 3]. There also exists a less well understood 
class of membrane invaginations, known as caveolae, that are less morphologically distinct 
than clathrin coated pits and resemble aspherical invaginations with a typical size of the 
order of 100nm[4]. There is good evidence that caveolae are involved in endocytosis [5] and 
play an important role in cell signaling[6]. Intriguingly they are probably also involved in the 
sensing of shear-stress [7, 8]. Such stresses would be expected to act on the tension 7 of the 
cell membrane. Thus, the present work may yield predictions for some of the direct physical 
consequences of surface tension, e.g. on inter-protein forces. In section II we discuss some 
qualitative effects of surface tension although a quantitative analysis of the effect of shear 
on membrane tension is beyond the scope of the present work. Recent elegant experiments 
demonstrate that better control of the surface tension may be achieved via micropipette 
aspiration facilitating studies of the effects of tension on membrane elasticity [9-11] and 
permeability to water [12]. 

It has now been shown that caveolins[13], a recently discovered class of membrane-bound 
proteins, are necessary for the formation of caveolae. These caveolins have a short membrane 
spanning sequence and N- and C-terminus polypeptide "tails", totalling about 150 amino 
acids, both found on the cytoplasmic side of the cell membrane. It is thought that these 
caveolin molecules are typically found in small aggregates of size b (a few nm) containing of 
approximately 15 molecules [14]. A crude thermodynamic (entropic) bound on the bending 
moment would be > k B Tb, higher for larger aggregates. 

Many theoretical studies have sought to calculate the effect of adsorbed polymer on flex- 
ible, fluid membranes, see e.g. [15-19]. One reason for this approach is that the simplified 
models available from polymer physics allow for a more or less exact computation of the 
entropic pressure exerted by the polymer on the membrane. Our work represents a general- 
ization of this approach. An anchored polymer is but one example, or model, of a membrane 
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inclusions that exerts a force on the membrane. Any membrane protein that isn't perfectly 
symetrical across the membrane should exert forces, the distribution of which will be dic- 
tated by the precise protein configuration. Furthermore, it can be argued that the physical 
effect of integral membrane protein can also be described by a well chosen force distribution. 
Hence there is a need for a general theory on the effect of arbitrary force distributions on 
flexible membrane. Our premier interest are large peripheral membrane protein such as 
caveolin. Our theory is, however, general enough to be applied to any flexible fluid interface 
with fixed tension, containing embedded impurities. 

Much of the previous work on interaction between membrane-bound objects has focused 
on membranes with vanishing surface tension [15-18], although there have been exceptions 
[20]. This work exists within the context of an extensive literature on membrane medi- 
ated interactions on tensionless membranes, see e.g. [21-23], which focussed mainly on 
Casimir-like forces originating from membrane fluctuations. It was argued that these could 
be either attractive or repulsive depending on the temperature and certain details of the 
model. The interest in tensionless membranes has its origin in the fact that isolated, self- 
assembled lipid membranes should be tensionless at equilibrium. Importantly this is not 
the case for many cell membranes. Other studies include tension as a (shape dependent) 
parameter determined self-consistently for closed membrane surfaces[24]. A study of the 
behavior of membranes under constant surface tension may shed light on the physics occur- 
ring on biological membrane, which are not truly at equilibrium and hence bear substantial 
tensions [25, 26]. Cells commonly adjust their surface tension to a set value via a mechanism 
known as Surface- Area- Regulation [2 7]. Hence membrane phenomena over sufficiently long 
timescales effectively occur at constant surface tension. Previous studies of the effects of rod- 
like objects embedded in fluid membranes[20] and impurities at gas-liquid interfaces [28, 29] 
share certain similarities with this study. 

II. ANALYSIS OF MEMBRANE DEFORMATION AND MEMBRANE MEDI- 
ATED INTERACTIONS 

Our aim is to construct a theoretical model for the membrane-mediated interactions 
between proteins bound to cell membranes. To achieve this we will denote the normal 
membrane displacement from its average flat state as u(r), where r is a 2D vector in the plane 
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of the unperturbed membrane. Motivated by a desire to analyze small membrane deviations 
we employ an expansion of the free energy of the membrane in powers and gradients of u. 
The free energy of a deformed fluid membrane, without inclusions, is well approximated 
by[30] 

F memb = J d 2 r g(VVr)) 2 + |(V M (r)) 2 ] (1) 

Which represents the truncation of an arbitrary expansion at order u 2 (odd orders are 
excluded by symmetry) and at the second derivatives of u. For a planar membrane, with 
normal in the indirection, the gradient operator is V = x-^ + The only two parameters 
that are needed to describe the physics of the membrane are k, the elastic bending modulus 
(typically [31] 20k B T for biomembranes) and 7 the surface tension (estimated to be [25] 1CT 1 
to 10~ 2 pN/nm). 

To include the work done by the pressure (normal force per unit area) /(r) exerted by one 
or more inclusions we exploit the fact that, for small deviations, the work done is merely the 
surface integral of the product of pressure and distance. The displacement has no effect on 
the pressure field at this order. Thus the total free energy, including the effect of inclusions 
is 

F = Jd 2 r [^(VVr)) 2 + |(V«(r)) 2 - /(r)u(r)] (2) 

where the term involving /(r) is chosen to have a minus sign to ensure that the membrane 
displacement has the desired sign. At thermodynamic equilibrium the free energy F must 
reach its minimum value, at which point u(r) satisfies the Euler-Lagrange equation 

£ 1 £ 2M ( r ) _ M = (3) 

K 

with C\ = V 2 , £2 — V 2 — k 2 and k 2 = 7//C indicating an intrinsic length scale k~ l m 
30nm (with perhaps 6nm < k^ 1 < lOOnm). We solve this under the boundary condition 
linir^oo Vti(r) = and that lim r _ >0 w(r) exists. Since the commutator [£i,£ 2 ] = it can be 
shown that 

u(v) = I J [G^r, r') - G 2 (r, r')] f(r>)d 2 r>, (4) 

where G\ and G2 are the Green's functions corresponding to the self adjoint problems 
determined by C\ and £2 (and their boundary conditions) respectively. We now seek to find 
the Green's functions of these self adjoint operators. 
For Ci we solve 

AGi(r, r') = V 2 G 1 (r, r') = -5(r - r') (5) 



The only boundary conditions compatible with this equation are Neumann conditions at 
infinity i.e. linir-r'^oo VGi(r — r') = 0. This corresponds to the membrane being asymptot- 
ically flat and yields the well known Green's function of Poisson's equation in 2-D 

G 1 (r,r') = --LlogA|r-r'| (6) 

where A is an arbitrary constant with dimensions of (length) -1 . We choose A = k for 
convenience, this choice corresponds only to a definition of the zero of the free energy. 
For C 2 we must solve 

(V 2 -A; 2 )G 2 (r,r') = -5(r-r / ). (7) 

subject to the boundary conditions specified above. This self adjoint problem has linearly 
independent solutions K (k\r — r'|) and I (k\r — r'|). These are, respectively, the zeroth 
order modified Bessel functions of the first and second kind. The function Iq(x) diverges as 
x — > oo therefore we write the Green's function as 

G 2 {v-v')^^K,{k\v-v'\) (8) 

This gives us the solution of equation (3), subject to the stated boundary conditions as 

«(r) = J G(r,r')/(r')dV (9) 

with the Greens function 

G(t, r') = - J- [K (k\r - r'|) + log k\r - r'|] (10) 

Z7T7 

Equations 9 and 10 represent a solution for the equilibrium membrane displacement u(r) 
due to an arbitrary force distribution. The formal analysis is concluded in the appendix, 
where we show that this solution is unique. 



Interaction Potential between two proteins 

Consider a fluid membrane with two, not necessarily identical, inclusions bound to it. 
Denote the pressure fields acting on the membrane surface due to each inclusion by tp^ and 
■?/>( 2 ). Without loss of generality we may let ip^ be centered at the origin and let the vector 
r joining the centers of the pressure distributions lie along the x-axis. Mathematically this 
means that we write the total pressure distribution as 
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/(r')=^(r / )+^ (2) (r , -r), (11) 

where ip^(r' — r) is the pressure distribution of the second inclusion centered at r, rather 
that the origin. Hence the Free Energy at equilibrium is given by F = F^ + F^ + <3>(r), 
where F^ is the self energy of the i th inclusion and $(r) is the interaction energy. We are 
interested only in $(r) as this is what determines the physics of interest to us here. By 
direct analysis or analogy with electrostatics this may be expressed as 

$(r) = -J d 2 r' J rf 2 r'> (1) (r')^ (2) (r /, )G'(r-r / + r") (12) 

where G(r) is the real space Green's function. 

For the case that we have two inclusions with circular symmetry we are able to determine 
the interaction potential exactly. Consider two circularly symmetric inclusions i.e. ^w(r') — 
^W(|r'|) = V W ( r ') such that ^(r 1 ) = for r' > b, where b is some maximum radius of 
the inclusion. Also each pressure field is such that no average force acts (/ ^(r')dV = 0). 
For inclusions which are not anchored to external structures, such as the cytoskeleton, this 
condition must be satisfied on general grounds. Indeed Newton's third law requires that any 
average force applied to the membrane by an inclusion must be equal and opposite to the 
reaction force applied by the inclusion to some external structure. If there is no external 
structure for this reaction force to act against then there can be no average force on the 
membrane. We assume that the pressure fields do not overlap so that r > 2b. Under these 
assumptions the logarithmic term in Eq (10) gives a vanishing contribution to the integral 
Eq (12) due to Gauss' Theorem. It corresponds to the electric field due to two circularly 
symmetric charged discs in 2-D, each carrying no average charge and thus with zero field in 
the region outside both discs. This analogy will be revisited when we discuss non-circularly 
symmetric inclusions in section IV below. The potential between two circularly symmetric 
inclusions is therefore given by 

$(r) = — / r'dr' / r'^rWOV^V) / d&l d9"K (k\r - (r' - r")|), (13) 
2tt7 Jo Jo Jo Jo 

the upper limits of integration are both b since ifj(r') = for r > b. We make use of the 

identities 

r K (k\r - (r' - r")|)<J0" = 27rK (kr)I (k\r' - r"|) (14) 
Jo 
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valid for r > |r' — r"|, and 



I (k\r' - r"\)d6' = 2irI (kr')I (kr") 



(15) 



valid for r' 7^ r" . Thus we are able to derive an exact analytical form for the interaction 
potential between non-overlapping (r > 2b) circularly symmetric inclusions 



where = 2ir J °° r'ip^(r')I (kr')dr' characterizes the strength of the i th mem- 

brane/inclusion coupling. By employing the expansion of Io(kr') the quantity £W can be 
thought of as a series of moments of the force distribution. The radial force is, as usual given 
by the derivative of the potential f r = — <9$/<9r. The interaction is everywhere repulsive. 
Some previous studies have reported attractive interactions between polymer chains grafted 
on tensionless membranes [15, 16]. In these studies the position -u(O) of the chain grafting 
point was fixed, effectively by a Lagrange force which then acted to ensure that there was 
no average force, although this condition was not identified explicitly by these authors. In 
the case of vanishing tension it was found that there was no interaction between inclusions 
unless they "touch" , i.e. their force fields overlap. One can understand this as follows: For 
tensionless membranes the membrane deformation under the force field is dictated by the 
minimisation of an energy similar to Eq (2) but with 7 = 0. The membrane deformation 
outside the extent of the force field simply minimises the curvature energy, hence has a zero 
mean curvature and does not contribute to the total energy. As two inclusions approaches 
each other, the total energy does not change until one force field couples to the membrane 
deformation directly under the other force field, since the deformation energy of any piece of 
membrane that it not in direct contact with a force field is zero. There is no characteristic 
lengthscale k^ 1 to give the interactions finite range when 7 = 0. In consequence, inclusions 
interact only when the two force fields overlap, in the case of grafted ideal polymer chains, 
the energy is decreased when the force fields overlap because they both want to deform the 
membrane in the same direction, resulting in an attraction. We note that the interaction 
between overlapping membrane inclusions is very sensitive to the direct physical, rather than 
membrane mediated, interactions between them, unlike Eq (16) which is universal. In the 
7^0 limit Eq (16) gives vanishing interactions, entirely consitant with these earlier results 
for tensionless membranes. 



$(r) 
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As mentioned in the introduction recent experiments on the effect of shear forces on living 
cells provide a crude control on membrane tension. In future studies we plan to calculate 
the precise effect of this on, e.g. surface phase equilibria and budding phenomena. 

III. MICROPHASE SEPARATION AT THE SURFACE OF CAVEOLAR INVAGI- 
NATION 

Membrane-mediated interactions between membrane-bound inclusions have been exten- 
sively studied theoretically [21] with, however, relatively little discussion of how experi- 
ments could easily support these predictions. Although presumably long range, these in- 
teractions are expected to be quite small, and are probably often dominated by stronger, 
(bio) chemically specific, short range phenomena. Long range interactions can, however, pro- 
foundly influence the phase behaviour of certain membrane proteins even if stronger, short 
range interactions, are present. As we emphasise below studies of such surface phases might 
provide an indirect probe of the membrane mediated interactions of interest here. This 
might be of particular relevance for the study of the phase behaviour of the protein caveolin, 
which is found at high density on certain membrane invaginations called caveolae (see sec. I). 

One peculiar feature of the caveolae bulbs is their texture. Distinct striations are seen at 
the surface of these buds[4] (see also Fig. 13.48 in [1]). These are now thought to correspond 
with the organization and alignment of caveolin oligomers on the membrane [32]. The 
observation of these surface stripes is intruiging, and the reader may find it interesting that 
radially symmetrical oligomers can give rise to non-symmetrical surface phases. We discuss 
this below, arguing that the stripe phase may be a signature of the membrane mediated 
repulsion between protein aggregates, such as those calculated in section II. Indeed, studying 
the phase behaviour of membrane inclusions may provide one of the best ways to test 
theory against experiment in this field. Molecular dissection of the caveolin protein has 
shown that caveolin oligomers strongly attract each other through contacts of specific protein 
sequences [6] (a short range attraction). It is well known that a solution of particle interacting 
via hard core repulsion and a short range attraction undergo a gas/liquid phase separation. 
This results in large dense regions (the liquid) coexisting with less dense regions (the gas). 
The caveolin oligomer being anchored to the cell membrane, there exists an additional 
membrane-mediated, longer range repulsion between between them, as demonstrated in 
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section II, which allow for a more complex phase behavior [44]. It has been recently argued 
at the light of computer simulations [33-35] that a long range repulsion should break the 
gas and liquid phases into microdomains. The short range attraction still locally drives 
a gas/liquid phase separation, but large aggregates are costly because of the long range 
repulsion. The microdomains are circular liquid islands at low density, circular gas regions 
at high density, and stripes otherwise. 

To gain a quantitative understanding of how stripes may form on caveolae we study 
the stability of an homogenous (2-D) solution of oligomers of area s b = irb 2 interact- 
ing via a potential V(r), which include a short range (the particle size b) attractive ex- 
ponential interactions and the longer range membrane mediated repulsion of Eq (16): 
V(r) = -E a e' r l b + E r K (kr), where E a and E r are respectively the strength of the at- 
traction and the repulsion. We look at small perturbations 8(p (with the conservation rule 
J dS5<fi(r) = 0) around the average surface coverage <f> . The free energy includes the pair in- 
teraction and the translation entropy of the inclusions, for which we use the gas-on-a-lattice 
model. Expanding the free energy: 

F = J ^(01og0+(l-0)log(l-0)) + i J ^0(r)T/(|r-rW) 
^ fdS 5(f) 2 1 [dSdS' ± ,^ r ,. 

5F = / — (T 7, S + o / ~2 8<f>(r)V(\r - r'\)5<t>{r') (17) 
J Sb 20o (1 - 0o) 2 J s z h 

The Fourier transform 5<f>(r) = S J (f^2^0g e * qr allows us to investigate the formation of 
structures. It leads to 

1 ^ 2 f d 2 n T 1 

sf = t;— tt^MI 2 with V q = — + -V q (18) 



2s b J {2nf yi ^ q 0o(l"0o) s b 

If V q , the Fourier transform of V(r), is sufficiently negative (attractive) for a given mode q 
then V q < and the mode is unstable. Substituting the expression for V(r), we obtain 



V Q = 2k B T 



+ ~E a - . „ XOXQ/0 + E r 



(19) 



_20 o (l-0o) V a {l + {bq?f/ 2 r k 2 + q 2 / 

In the absence of long range repulsion, the most unstable mode is always q = 
(macrophase separation), and the liquid/gas transition is observed provided 2E a (p (l — O ) > 
k B T. Because of the membrane-mediated long range repulsion, the function V q presents a 
minimum if E r > ^E a (kb) 4 ~ 1CT 3 -E a (for the typical numbers b ~ 5nm and k" 1 ~ 30nm). 
One sees that although the strength of membrane-mediated interaction (Eq (16)) is expected 
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to be quite small: E r ~ lO~ 2 k B T (see sec. V.A, or ref.[36] for a more detailed analysis), it 
is of much longer range, and should be competent to produce a well ordered phase. 

The phase separation occurs preferentially for a mode given by dV q /dq = 0, and periodic 
arrays of dense and dilute regions are observed. In the limit kb ^ 1 the structure have a 
typical size 2%/q* which is independent of the range of the repulsion: q*b = {2E r /'iE a ) 1 ^ ~ 
\. It thus defines dense stripes of width about 5 particle diameters, which agrees quite well 
with the experimental observations. The existence of a stripe phase at the surface of the 
biological structures known as caveolae is thus quantitatively consistent with the existence 
of a membrane-mediated repulsion between proteins described by Eq (16). 

This analysis does not help us to discriminate between circular domains and stripes. 
Computer simulation have clearly shown that this transition indeed exists for high enough 
surface coverage [34, 35]. 

Note that the formation of striped mesostructures in membrane have been predicted on 
the basis of curvature effects only[37, 38]. However, it can be shown that their existence at 
the surface of invagination of such large curvature as the caveolae (R ~ 50nm) is inconsistent 
with the existence of invaginations with a well defined size. 



IV. FAR FIELD INTERACTIONS OF NON CIRCULARLY SYMMETRIC IN- 
CLUSIONS 

If the force distribution of the inclusions is not circularly symmetric then no general 
analytic solution to Eq (12) exists. However we can still proceed by examining the far-field 
interactions of slightly asymmetric inclusions. In order to do this we will parameterize the 
asymmetry of the s th inclusion by 

£>« = J dW/^\r') = ^ " + e(S) c j (20) 

where i and j are cartesian indices, c controls the average (isotropic) magnitude and e 
controls the anisotropy of a slightly "elliptical" force distribution, extended in the x-direction 
and contracted in the ^/-direction (for e > 0). This force distribution can be rotated by 
the usual rotation matrix R(9) to obtain the effective moment Ri m (6^)Rj n (9^)Dmn of an 
"elliptical" force distribution rotated by an angle 9^ relative to the x-axis. This expression 
involves two rotations, one for each factor of r' in Eq (20). 
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It is important to first note is that the interaction can be separated into two terms, one 
that is an integral involving K (k\r — r' + r"|) and the other that is an integral involving 
log k\r — r' + r"|. The first of these terms gives a contribution to the interaction potential that 
is dominated by a term like Eq (16), with corrections due to the asymmetry that are smaller 
(and also asymptotically exponentially short ranged). This conclusion is independent of 
the precise choice of asymmetry, provided only that it is small. The contribution from the 
second term is more interesting and can be shown to give rise to a small but long ranged 
correction 5$ to the interaction potential. This can be best understood by way of multipole- 
like expansion of log k\r — r| where |f| <C |r] and we choose the direction of r = rx to define 
the x-axis (the 1 direction), without loss of generality. Thus f\ and f 2 are the x and y 
components of f and 

log k\r - ^ log kr - f -l + ^1 - fLfffl + 6f m "J? - f i + o (1) (21, 

The conditions of no overall force J d 2 r'^(r') = and no overall moment J d 2 r'r'^(r') = 
for each inclusion mean that the first non- vanishing contribution to 5$ ~ l/?" 4 , the classical 
result for electrostatic quadrupole-quadrupole interactions in two dimensions. Using f = 
r' — r" in Eq (21) to leading order in 1/r we have[45] 



5$ = _1_ J d 2 r' J dVV ( V)^ (2 V)log^|r-r' + r"| 

= I 6 \ cos2(gW + ^ 2 )) (22) 

This result demonstrates that elliptical inclusions attract if arranged so that 0^ + 0^ = nit 
with n an integer, as reported in an earlier study of rod-like inclusions [20]. This condition 
represents a degenerate family of orientations in which the orientation of the quadrupoles 
has reflectional symmetry about the mid plane. This perturbative interaction is long-ranged, 
scaling like 1/r 4 . Interestingly this dominates the exponentially short ranged repulsion from 
Eq (16) for large enough separations. 

There is a straightforward way to understand the appearance of an algebraic potential 
for anisotropic inclusions. The physics of Eq (5) and is that of the Poisson equation in 
which the force distribution is analogous to the electrostatic charge. The potential due to 
a point charge in 2D is therefore Eq (6). Gauss law tells us that there is no field outside 
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a circularly symmetric charge distribution that has no overall charge (imbalance). This 
condition is analogous to our requirement here that there be no overall force J d 2 r'f(r') = 
by Newton's law. If there is no field the potential is constant (zero without loss of generality) 
and no forces act on the inclusions. However, if we relax the condition that the force (charge) 
distribution is circularly symmetric there is no longer a simple symmetry argument that the 
radial and azimuthal components of the field must be zero. Indeed they are not. As is usual 
in multipole expansions in electrostatics an algebraic potential results. It is here analogous 
to the quadrupole-quadrupole interaction in electrostatics, since there are no dipole moments 
J d 2 r'r' /(r') = 0. This is due to the fact that there can be no external first force moment 
(torque) on the membrane if it is not anchored to any external structure against which the 
moment can act. 

The physical consequences of these interactions are potentially significant. In spite of 
the fact that the dominant interaction is repulsive the surface concentration of membrane 
inclusions may often be maintained at fairly high surface fractions by the regulatory mech- 
anisms of the living cell. The composition of the surface coats of caveolae may be further 
enriched in several important membrane components, including cholesterol and proteins in 
the caveolin family, as well as others[6, 13]. This means that there is an effective surface 
pressure, driven by a chemical potential difference that forces the inclusions to partially 
overcome their repulsion. In this situation the attractive 0(e^e^) interactions might be- 
come significant, leading to an in plane anisotropic phase separation such as is observed in 
model systems[33-35, 39] and on the coats of caveolae[4]. For caveolae the typical inter- 
actions seem too weak [36] to be solely responsible for the stripe morphology, which may 
rather be dominated by specific attractions, see section III [32]. Finally the existence of 
an attractive interaction that can arise from fluctuations in the inclusion force distribution 
(shape), as parameterized by e, suggests the possibility of an attractive, fluctuation driven 
force reminiscent of Van der Waals forces. This mechanism should extend to systems such 
as wetting droplets at gas-liquid interfaces. These can potentially exert much stronger forces 
on the interface, resulting in a more exaggerated effect (see section VB). 
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V. ORIGIN OF THE FORCE DISTRIBUTION 



Up to this point we have merely postulated the existence of a force distribution ip(r') due 
to a membrane bound inclusion. This force can arise from direct mechanical effects, such as 
due to the geometry and shape of the inclusion, e.g. the protein clathrin, or entropic effects 
due to the asymmetric anchoring of, e.g. flexible hydrophilic polymers, to a membrane. The 
budding of caveolae is now thought to be driven by the protein caveolin which is known to 
form small oligomers resembling a number of polypeptide chains extending only from the 
cytoplasmic side of the membrane[14, 40]. 

We have in mind that we have demonstrated that interactions should exist in general, 
but haven't yet adressed the question of how big they might typically be. We are particu- 
larly interested in any models that can be chosen to be broadly comparable with membrane 
proteins. In this section we will present two examples of models for the origin of a force 
distribution arising from anchored polymer chains. This approach represents the applica- 
tion of a simple, and hence rather "idealised", theory drawn from polymer physics. Its 
exact quantitative applicability to polypeptide chains is almost certainly limited to a small 
subset of peptide sequences which either resemble a diffuse random coil or a dynamic, dense 
hydrophobic globule. Nontheless in the absence of any detailed microscopic information 
on a specefic inclusion's configuration and stability the approaches that we outline below 
probably represent the best chance for us to obtain a rough idea of the likely scale of the 
effects that we have described in this article. It is quite possible that the forces exerted by 
membrane proteins with significant well defined tertiary structure may be somewhat higher 
than the estimates that we will present below. This encourages us to consider our estimates 
as approximate lower bounds on the scale that these effects might reach. 

In polymer physics terminology the two models that we consider for the force applied 
by a membrane-bound polymer will correspond to the two limiting of "good" and "poor" 
solvent conditions respectively [41]. We will discuss below how these models might apply to 
caveolin homo-oligomers. Other recent studies of the physics of flexible polymers grafted 
onto tensionless fluid membranes [15-1 8] undoubtedly share similar motivation. One feature 
that these earlier studies all have in common is that there is only an attraction for polymer 
force distributions that overlap one-another on the surface of the membrane and, when this 
is the case, the interaction is attractive. This is opposite in sign to the interactions that 
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we predict on membranes that bear tension, see Eq (16). Furthermore our interactions are 
extended, with ranges ~ A; -1 . 

A. Anchored polymers and polymer aggregates: Good solvent conditions 

Motivated by possible biophysical relevance we investigate the effect of flexible polymer 
chains anchored to a small patch of membrane, of radius a. We assume that the chains are 
in a good solvent, i.e. that they are found in extended, hydrated random coil configurations. 
The chains form a hemisphere of outer radius b > a, see figure 1. This model may be used to 
treat polymer homo-oligomers made up of a general number Q of chains. As such it shares 
many features with those that are known for caveolin homo-oligomers[14]. 

Within the "corona" of the polymer hemisphere, at radial distances b < r < a from the 
center of the polymer aggregate, there exists a characteristic correlation length £(r) for the 
polymer chains which crudely, represents the spatial distance between collisions between the 
segments of a chain with other chains or the membrane [41]. The conservation of number 
of chains implies [42] ^Anr 2 ~ Qi{r) 2 for a < r < b. This is since the surface area of a 
hemisphere of radius r is filled by close packed blobs, up to a constant of order unity. From 
this one may immediately deduce the scaling of the correlation length £(r) = 2 r - The 
work done in generating each blob is k B T, independent of the blob size. Thus we may write 
the pressure in this region as the energy per blob divided by the volume of a blob [42] 

k B T (Q\^k B T 

/(r) = Wf = w (23) 

for a < r < b. For r > b the pressure is here assumed to be zero [46]. The pressure in the core 
binding region is assumed to be constant and must involve a total force equal and opposite 
to that applied by the corona so that no average force acts, as required by Newton's second 
law. 

E a~ 3 if < r < a 



f(r) = | if a<r<b (24) 

if r > b 

where E Q = f(0)a 3 = 2k B T(l — a/b)(Q /2ir) 3 / 2 is a characteristic energy related to the force 
applied at the center of the aggregate. This model represents one approximate framework for 
the understanding of forces applied by membrane biopolymers, see Fig 1. It possesses several 
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important features: It involves no average force being applied to the membrane. It involves 
no average first force moment by symmetry. The first non-vanishing moment is the second, 
corresponding to a bending moment with magnitude controlled by f r 2 f(r)d 2 r ~ E Q a with 
the characteristic energy scale E Q > k B T that is entropic in origin. 

Models similar to this have been proposed elsewhere[15-18] although, as emphasized 
previously, the physical behavior of these inclusions on fluid membranes that are under 
tension is rather different. 

We now proceed to give an estimate of the membrane deformation and interaction po- 
tentials due to membrane-bound biopolymers within the good solvent polymer brush model 
described above. 

For biological membranes primarily comprised of a bilayer of phospholipids the tension 
and rigidites are typically[25, 31] 

k = lOk B T - A0k B T k Uo = 20k B T 7 = 1(T 5 - 1(T 3 J/m 2 lhio = 10~ 4 J/m 2 (25) 

where xuo refers to what we would take as a single "typical" value for biological membranes. 
From this we obtain the lengthscale 

k^ 1 = 6nm — lOOnm k^ o = 30nm (26) 

which may crudely be identified with the range of the interactions. 

1. Membrane deformation due to the inclusion 
The deformation at r = is 

Eq 

u b rush(r = 0) = — a (27) 

K 

in a frame where u = at infinity. This therefore represents the total magnitude of the 
normal deviation of the membrane expressed in terms of a dimensionless ratio and a 
characteristic length a. 

For parameters that might be typical of caveolin homo-oligomers, a = 2nm and b = 5nm, 
we find 

a(k~ l = 6nm) = 0.7nm a(k~ l = 30nm) = 1.5nm a(k~ l = lOOnm) = 2nm (28) 
much the same as for slightly larger inclusions with a = Inm and b = lOnm 

a{k~ l = 6nm) = 1.3nm a^k' 1 = 30nm) = 2.7nm a{k^ = lOOnm) = 4nm (29) 
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2. Membrane mediated interaction energy between inclusions 



The interaction energy is expressed as 

^brush = ^P (30) 
K 

which is the product of an energy Eq/k, typically of the order of k B T, and a dimensionless 
number (3. For a = 2nm and b = 5nm 

Pik' 1 = 6nm) = 3 KT 2 P(k~ x = 30nm) = KT 3 p(k^ = lOOnm) = KT 4 (31) 

while for a = lnm and b = lOnm we find instead 

Pik- 1 = 6nm) = 0.13 p(k~ l = 30nm) = 5 10~ 3 P{k~ l = lOOnm) = 5 10" 4 (32) 

B. Poor solvent conditions 

If the polymer chains are more hydrophobic, and the solvent conditions are poor, the 
polymer chains may collapse into a tight, roughly hemispherical, region from which water 
is largely exuded. This represents the natural opposite limit to the good solvent polymer 
model considered above. In the poor solvent case we propose to model the inclusion as 
if it were a fluid droplet that partially 'wets' the membrane surface[43]. At any point on 
the circumference of the drop the forces per unit length (due to the surface energy) are 
in equilibrium. Resolving these force perpendicular to the membrane we find a total force 
acting upward 

F = 2nasm9 (33) 

where b is the radius of the drop, a is the energy cost per unit area of producing the interface 
between the inclusion and the external poor solvent, typically water, and 9 is the contact 
angle given by Young's law [43]. Since we are in equilibrium this force must be balanced by 
the force acting on the interface between the drop (inclusion) and the membrane for r < b. 
To find the corresponding pressure we need only consider the partial wetting droplet as if it 
formed part of a larger sphere. The pressure is constant everywhere inside the sphere and 
is given by the Laplace law 

„ BE 2a 

P = — — = — , 34 
dV R' K J 
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where R — b/ sin9 is the radius of the sphere. The radius R can therefore be related to 
the volume of the drop (and b) by simple geometry. Hence we may write the pressure 
distribution as 



where 6(x) is the Heaviside unit step function. 

This model represents a different possible physical origin for the pressure distribution that 
might be valid, e.g. for surface-anchored hydrophobic polymers in the poor solvent regime, 
when little water penetrates the protein chains, see Fig 2. It also preserves the same features 
as the good solvent polymer model: It involves no average applied force or first force moment 
and has a bending moment f r 2 f{r)d 2 r ~ E Q b that is enhanced by the localization of part of 
the force distribution around the exterior circumference r = b of the inclusion's footprint on 
the membrane. The characteristic energy scale is now chemical in origin, as it is controlled 
by chemical parameters such as a. It is difficult to give a quantitative scale for this although 
typical oil water interfacial tensions of o ow > 3 x 10~ 2 Jm -2 suggest E > u ow b 2 180k B T 
with b ~ 5nm. This large energy suggests thermodynamically "strong" interactions. 

1. Membrane deformation due to the inclusion 

For b = lOnm, a = AQdyn / cm[A7] the membrane's normal deformation at r = is 

Uo(k~ l = 9nm) = 37nm Uo{k~ 1 = 30nm) = 65nm Uo(k~ l = 90nm) = 92nm (36) 
for b = lOOnm a = 40dyn/cm 

Uo(k^ x = 9nm) = 2700nm Uoik" 1 = 30nm) = 15/im Uo{k~ 1 = 90nm) = 37/im (37) 

2. Membrane mediated interaction energy between inclusions 
The interaction potential energy has a simple analytical form: 



/(r) = -|(6(J(r-6)-26(6-r)) 



(35) 



$ = 




(kb) 2 K (kr)(3 2 



(38) 



2-KK 



with 




(39) 
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For b = lOnm and a = AOdyn/cm the energy scale E (simply the potential $ divided by 
K {kr)k B T) is 

Eoik- 1 = 9nm) = 25 E Q (k~ l = 30nm) = 2 E Q (k~ l = 90nm) = 0.2 (40) 
for b = lOOnm a = 40dyn/cm 

Eoik- 1 = 9nm) = 10 9 E Q (k~ l = 30nm) = 3 10 6 E^k' 1 = 90nm) = 2 10 5 (41) 

VI. CONCLUSIONS 

We show that a model for the interactions between inclusions bound to fluid membranes 
can be solved exactly for circularly symmetric inclusions bound to membranes that are under 
tension. We argue that proteins bound to cell membranes can apply no average force, or first 
moment of force, to the membrane unless they are also anchored to an external structure, 
such as the cytoskeleton. By this we mean that the total of any and all forces that the 
membrabe protein acts to apply downwards against the membrane must exactly equal the 
magnitude of similar forces applied upwards. This is a consequence of Newton's laws of 
motion. 

We proposed idealized models for the origin of these forces due to either entropic or 
direct chemical interactions. The interactions between two circularly symmetric inclusions 
are repulsive and are asymptotically exponentially short ranged with a typical extent on 
biological membranes that is of the order of A;" 1 ~ 30nm. This result contrasts with the 
attraction predicted to appear between similar inclusions on tensionless membranes. For non 
circularly symmetric inclusions we predict an additional algebraically long-ranged attractive 
contribution. 

We discuss how competing attractive and repulsive interactions are known to sometimes 
produce a stripe-like morphology. These phases are reminiscent of the stripes observed on 
the surface of caveloae, which we refer to as "gnarly buds", and we discuss how a similar 
lateral phase separation could occur in these systems with typical length scales comparable 
to the ranges of the interactions. 

We believe that the results discussed here may have wider applications in understand- 
ing biological phenomena, including lateral phase separation in phospholipid membranes 
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and endocytotic budding. Indeed, a detailed theoretical analysis of the formation of cave- 
olae invaginations on model phospholipid membranes utilizes many of the results derived 
here [36]. 
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Appendix: Uniqueness of the solution for the equilibrium membrane displacement 
due to an arbitrary force distribution 

Equations 9 and 10 permit the calculation of the equilibrium membrane displacement 
due to an arbitrary force distribution. This solution is unique, as may be demonstrated as 
follows. 

We may add to our solution any solution of the homogeneous version of equation (3) 
without changing the result. Let ui(r) satisfy 



and the boundary conditions imposed. Multiplying by Mi(r) and integrating we have 



Integrating by parts and noting that the integrated terms are zero by virtue of the boundary 
conditions we obtain 



The integrand is everywhere positive so we must have V 2 wi(r) = and Vtti(r) = every- 
where. Thus the only admissible solution of the homogeneous equation is u(r) = constant 
this correspond only to a redefinition of the zero of the displacement and so without loss 



V 2 (V 2 -A; 2 ) Ml (r) = 



(42) 




(43) 




(44) 
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of generality we take this constant to be zero. This freedom of choice is a consequence 
of the translational symmetry of the problem in the direction perpendicular to the planar 
membrane. 
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approaches the oil-water tension a = Wdyn/cm — AOdyn/cm = 2.5 — lOfesT/nm 2 , and a 
curvature radius for the surface of the drop R = 1 — lOOnm. These combines to give a pressure 
fb = cr/R = 2.5 10~ 2 — 10/csT/nm 3 . We arbitrarily choose b = R/2 to give an idea of the 
scales and define the energy unit E b = f b b 3 = aR 3 /8, so E b = 0.3k B T - 10 4 k B T. With 
&bio = ^Odyn/cm = lO/cB^/nm 2 and Ru = lOnm we obtain f b ,bio = lfesT/nm 3 and E b)b i = 
125/cbT. We give only the results for a = AOdyn/cm and for k = 20/cbT for compactness. 
With thses values /c -1 = 9 — 90nm. 
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FIG. 1: Sketch of the form of the anchored polymer aggregate in good solvent conditions, when 
the polymer configurations are somewhat extended from the membrane. Our model for the force 
distribution f(r') is given by Eq (24) and represents the force per area applied to the membrane 
by the polymers. The membrane is pushed down by the "corona" of the grafted polymers out to 
r = b = 5nm and is pulled upwards by the anchored "core" inside r' = a = 2nm as a result. For 
aggregates residing on the cytoplasmic face of the membrane, including caveolin homo-oligomers 
the cell interior would be above the membrane. The resulting deformation u(r') is shown out to 
r' = 15nm for the following values of the parameter values; 7 = lCT 4 J/m~ 2 , k = 20k B T and 
Q = 15. The force and deformation can be expressed in terms of the characteristic energy E , see 
Eq (24). 
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FIG. 2: Sketch of the shape of the anchored aggregate in poor solvent conditions, when the polymer 
forms a droplet on the membrane that is assumed to largely exclude water. Our model for the force 
distribution f(r') is given by Eq (35). The membrane is now pulled upwards by the resolved Young's 
force on the contact line at r' = b = 5nm and is pushed downwards by the Laplace hydrostatic 
pressure for r' < b as a result. The resulting deformation u(r') is shown out to r' = 15nm for the 
same total integrated force J d 2 r'\f(r r )\ shown in Fig 1, for comparison. All the other parameters 
values are also as given in the caption to Fig 1. Even for the same integrated absolute force the 
membrane deformation is an order of magnitude larger than shown in Fig 1, as well as having 
opposite sign. The enhanced effect is due to the concentration of the forces at the exterior of the 
aggregate. 
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